Thermo-magneto coupling in a dipole plasma 
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On a dipole plasma, we observe the generation of magnetic moment, as the movement of the levitating 
magnet-plasma compound, in response to electron-cyclotron heating and the increase of /? (magnetically- 
confined thermal energy). We formulate a thermodynamic model with interpreting heating as injection of 
microscopic magnetic moment; the corresponding chemical potential is the ambient magnetic field. 

PACS numbers: 



The RT- 1 device confines a high-temperature (electron tem- 
perature Te ^ 10 keV) plasma in a dipole magnetic field that 
is generated by a levitating superconducting magnet Dll-H] ; see 
Fig.[T] When a high-beta (local /3 ^ 0.7) plasma is produced, 
we observe an appreciable amplitude of vertical motion of the 
levitating magnet-plasma compound, while the magnet posi- 
tion is regulated by a feedback control system [5]. Interpreting 
this phenomenon form thermodynamic view point, we will de- 
lineate an interesting property of magnetized plasmas. 

Let us start by analyzing mechanics. We denote by z the 
vertical displacement of the magnet from the equilibrium po- 
sition (we use r-O-z cylindrical coordinates). From the time- 
series data of the coil position and the controlled current (/l) 
in the lifting coil, we can estimate the change of forces acting 
on the levitating magnet-plasma compound: denoting by M 
(=112 kg) the mass of the magnet (the mass of the plasma is 
ignorable), the equation of motion in the vertical direction can 
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be written as 






(1) 



On the right-hand side, Mg is the gravity and Fm is the mag- 
netic force in the vertical direction: 

F,n := f e,-{Jx Bl) d?x « -2^RIBL,r. (2) 

where Bl is the magnetic field applied by the lifting mag- 
net {BL,r is its radial component), J is the current density 
in the magnet-plasma compound, and / is the total current 
in 6'-direction. Here we have approximated J by a ring cur- 
rent of radius R. Invoking the conventional magnetic moment 
A^ = ttR^I, we may write (using V ■ B = r^^d{rBr)/dr + 
dB^/dz = and dBr/dr « Br/r) 
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We define M to be positive and B^r to be negative {Blz and 
dB^z/dz to be positive), and then, F^ is positive (upward). 
We denote 



dBj 



dz 



(4) 



to write Fm = M.G- At the equilibrium point (z = 0), we 
define AA — AAq and Q = Go- The equilibrium condition 
reads as MqQq = Mg. 

While (O is derived for the conventional magnetic moment 
of a loop current, we may use it to "define" the total mag- 
netic moment of the magnet-plasma system. In what follows, 
we evaluate Q at the barycenter of the levitating magnet, and 
define J\A :— F^/G by the total magnetic force F,„ on the 
levitating magnet-plasma compound. Linearizing ([U in the 
neighborhood of the equilibrium point {z — Q, G ~ Go and 
M = Mo), we obtain 



.rd'^ ,, dG 



Mo-^dlL+GodM, (5) 
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FIG. 1: Schematic drawing of the RT-1 device. A dipole mag- 
netic field is produced by the levitating superconducting magnet. 
The field strength in the confinement region varies from 0.5 T to 
0.01 T. Plasma is produced and heated by ECH (8.25 GHz, 25 kW, 
and 2.45GHz, 20 kW systems). 



where {dG/dlL)dlL represents the variation of G due to a 
perturbation d/i in the lifting magnet (we define the sign of 
II so that d/i > increases Bl^ and G) [? ]■ The inertial 
force on the left-hand-side of Q can be estimated by the time- 
series data of the coil position. Evaluating the first and second 
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FIG. 2: Typical waveforms of (a) the vertical plasma position z 
(measured by laser position sensors), (b) the lifting-magnet current 
II (feedback controlled in response to the position signal), (c) (1) 
the inertial force, (2) magnetic force by perturbed d/L, and (3) the 
remaining term of (|5) corresponding to the magnetic force by the 
perturbed magnetic moment dA^, (d) the perturbed magnetic mo- 
ment normalized by the magnetic moment A^o of the superconduct- 
ing magnet, and (e) the normalized plasma energy (/? averaged over 
the plasma volume) estimated by diamagnetic signals. ECH is in- 
jected for 0.5 < f < 1.5 s. 



terms on the right-hand side of (|5) by measured dz (= z) and 
d/L, we obtain the remaining third term, by which we can 
derive AM. . 

We observe that the magnetic moment M. increases as the 
plasma is heated; in Fig.|2]we compare the waveforms of the 
change in hA and the volume-average /3 estimated by dia- 
magnetic signals (for the detail of the measurement, see |3l)- 
We may explain the increased magnetic moment in terms of 
the diamagnetic current driven by the plasma pressure. The 
unique structure of this device — a "levitating" confinement 
system — therefore provides us with a particular method of 
estimating the plasma pressure by measuring the mechanical 
motion of the magnet; in Fig. [3] we show an experimental re- 
lation between /3 and dA^ . 
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FIG. 3: Experimental relation between the plasma- volume-averaged 
/?avc (estimated by the diamagnetic signal) and the increment of the 
magnetic moment dA^ (normalized by the magnetic moment Alo of 
the superconducting magnet). The experimental scaling (dotted line) 
shows AM I Ma = 0.75/3ave 



The main theme of this brief communication, however, is 
not the practical application of the magnetic moment for mea- 
surement. Examining this phenomenon from a thermody- 
namic view point, we notice an interesting implication, and 
that is the subject of present practice. Injecting electron cy- 
clotron heating (ECH) power, we increase the internal (ther- 
mal) energy of the plasma (electrons). In the language of ther- 
modynamics, giving a heat 8Q causes a change AU of the 
internal energy U (here we denote a general variation by 6X, 
while a variation of a state variable Y is written as dF ); the en- 
ergy U is the combination of the thermal, mechanical, gravita- 
tional and electromagnetic energies, and 8Q, in general, may 
cause variations in every component of the energy, resulting 
in changes in macroscopic quantities including those mechan- 
ical (vertical velocity), gravitational (vertical position), and 
electromagnetic (magnetic moment). Writing the first law as 



AU = bW ^^ 



(6) 



the term SW represents whole such contributions from macro- 
scopic quantities to the energy balance. In textbook thermody- 
namics, we often assume that 6W = —PAV with a pressure 
P and volume V, and then, the coupling of the thermody- 
namic energy and the macroscopic mechanical energy is only 
through compressible motion of fluid. Needless to say, possi- 
ble processes are much more rich in a plasma. 

As mentioned above, we observe that heating SQ causes a 
change in the magnetic moment M and subsequent changes in 
the vertical position and (feedback controlled) lifting-magnet 
current. To describe the "thermodynamics" of this system, we 
have to formulate the relations among 5Q, AM, AIl, and Az. 
Here we proffer a "grand-canonical model" to understand this 
thermo-magneto coupling. We do not intend to challenge the 
aforementioned elementary understanding in terms of the dia- 
magnetic current. Instead, our new perspective will delineate 
an interesting property of a magnetized plasma in a more suc- 
cinct picture. 

The energy of a magnetic moment la] is at the core of the 
first law connecting the plasma, the magnet, the heating sys- 
tem, and the lifting system. When an external magnetic field 



Bl is applied, the magnetic moment Ai has a mechanical 
potential energy f^ 



K„ = -MBl. 



(7) 



where Bl is the average of ez ■ Bl over the levitating magnet- 
plasma compound. As well known [8J, the "total energy" of a 
magnetic moment, including the electric energies of the levi- 
tating and lifting currents, is —V„i, but we must use Vm to de- 
rive mechanical forces and corresponding works. Combining 
Vm with the gravitational energy Ug = Mgz and the kinetic 
energy Uk = p1/{2M) (pz = Mdz/dt is the momentum), 
we obtain a Hamiltonian 



H := Uk 



Ug + K, 



(8) 



The corresponding Hamilton's equation of motion reproduces 
(|5]l. The explicit dependences of Vm on the parameters I^ and 
TM yield changes of H: 



dH_ _ dH_ _ dVrndh dV^dM 
'df ^~dr~ diL dt ~^^M~dr' 



(9) 



As already remarked, Vm is not the right energy to be 
inserted into the first law; we have to add the electric en- 
ergies in the levitating and lifting systems, which amounts 
2A4Bl = -^2Vm 1 8]. Hence, the magnetic moment acquires 
an energy, when put in an external magnetic field Bl, 



Urn - Vm + 2MBl = MBl = -Vm. 



(10) 



Notice the flip of the sign of energy. In addition to this mutual 
energy, the magnetic field of the total system (consisting of 
Bd that is produced by the dipole magnet-plasma compound 
and Bl that is produced by the lifting magnet) has also the 
self-energy Us that may be written as Ja] 
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(11) 



where Bd is an average of e^ • Bl>- 

Including the thermal energy Ut of the plasma, the total 
energy of the system is 

U = Uk + Ug + Urn + Us + Ut. (12) 

The first law, combined with the "mechanical law" (|9]l, reads 



dU = MdBL + BLdM + dC/, + dUt 



(13) 



We find that MdBL = M{dBL/dlL)dlL contributes to the 
mechanical work dH (on the magnet-plasma subsystem) by 
{dVm/dlL)dlL = -[M{d'^BL/dlLdz)dlL]dz (notice the 
flip of the sign). On the other hand, dM is "caused" by heat- 
ing 5Q, thus we may relate the term BLdM with 5Q (the 
latter also includes energy loss). 

To delineate the relation between dA^ and 5Q, we invoke 
the microscopic magnetic moment ji — (meW^)/(2i?), where 



B is the local magnetic field in the plasma region, me is the 
mass of an electron, and Vc is the velocity of cyclotron motion. 
The power of ECH, first of all, increases v1 (and then, excites 
macroscopic processes). The perpendicular thermal energy 
Ut.i_ is the sum of Bjjlj over all particles (labeled by j = 
1, 2, • • • ). With an average magnetic field B, we write 



Ut. 



J2BH-13J2 



H- 



(14) 



In view of ( fT4l i. we may rephrase "heating" as injection of 
microscopic magnetic moments jij , and then, B is an effective 
chemical potential. 

To relate the microscopic magnetic moments /ij with the 
macroscopic one Mp (we denote by Mp the plasma's contri- 
bution to M), we put 



Mp^vY, 



fij 



(15) 



with a geometric factor T), which we can estimate as follows. 
By the levitating magnets's current /q and the length scale £ 
of poloidal magnetic field lines {£ ^ 27ra with a minor radius 
a), we estimate B = jiqIq/I. Normalizing by J\Aq = ttR^Iq, 
we obtain 



M, 



Ut. 



Mo TTR^m^/fio 



= v 



V 



2t:R^£ 



/3-L, 



(16) 



where V is the volume of the plasma and /3± :— 
Ut.±/{VB'^ /2fj,Q) is the average beta ratio of the perpen- 
dicular plasma pressure. On the other hand, we estimate 
Mp = irR^Ip, where Ip is the diamagnetic current induced 



by the perpendicular pressure P± 
we obtain 

M, 



Estimating Ip = £P±/B, 



Mo 



2 ' 



(17) 



FigureU shows a reasonable agreement. Comparing (IT&i and 
([TtI i. we estimate V = ■nR^£/V ^ R/a. Since the change of 
the superconductor's current in response to dz, d/^ or dMp 
is of second order, we may assume dA^ = dMp. 

Now we have a more explicit representation of the thermo- 
magneto coupling processes included in the first law ( fTJt : 
denoting by C/j y the remaining parallel component of the 
thermal energy Ut and by A^l the coefficient such that 
fBld^x/^IQ = MLBL, 
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Ml\ dBL 
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Bd B 
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dir. 
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t,\\ 



(18) 



The first term on the right-hand side (induced by d/i) is the 
process connected to the lifting magnet system. The second 
term (induced by dA^) is the "ECH heating" SQech (or, in 
our language, injection of magnetic moments); the component 
{B/'D)dM goes to the thermal energy Ut,±, while the other 



components change macroscopic magnetic energies Um and 
Us, as well as mechanical energies Uk and Ug (through the 
mechanical potential energy Vm — —Um), which we observe 
as the change of z. The remaining abstract terms dC/^ || (paral- 
lel energy change) and 5Q' (heat processes including thermal 
conduction, energy loss with particle transport, etc.) are not 
the direct subject of the present analysis. 

We have made an attempt to understand and interpret the 
observed macroscopic thermo-magneto coupling in a dipole 
plasma produced on the RT-1 magnetospheric device. The 
most abstract thermodynamic first law (|6) has been given a 
more concrete and dissected form (ITS^ that elucidates the in- 
ternal and external thermo-magneto processes; the conven- 
tional expression of ECH as heating 5Q has been rewritten 
an injection of magnetic moment clM., and its partition into 
different terms of energy has been specified. 

What is rather nontrivial is that a magnetic moment m is 
an axial vector (or, a pseudo-vector) having an odd parity; the 
M. is the z-component of m (i.e. m = Me^), which can be 
regarded as a pseudo-scalar. Multiplying m (M) by the other 
axial vector B (a pseudo-scalar B), we obtain a scalar that 
can be related to an energy or some thermodynamic potential. 
Remember that the enthalpy U + PV of a neutral fluid cou- 
ples with a product VP • u of two vectors VP and u (fluid 
velocity). Or, more simply, we write the work as PdV (or 
—Vc\P for estimating enthalpy) with two scalars P and V . 
Relating the pressure P to the thermal energy by an equation 
of state, we can close a thermodynamic relation. To describe 
a thermodynamic model of a plasma, therefore, we have to 
find a relation between an axial vector (pseudo-scalar) and the 



thermal energy — there must be an intrinsic mirror-symmetry 
breaking to make such a relation possible. We have proposed 
a "grand-canonical model" with a pseudo-sealer chemical po- 
tential (that is the ambient magnetic field introducing the sym- 
metry breaking). 

We end this brief communication with a comment to ex- 
tend the scope of the paradigm of pseudo-scalar chemical po- 
tentials; different mechanisms of magnetic field (axial vector) 
generation can be related on a unified perspective. Remem- 
ber that the helicity K :— / A • Bd?x is also a pseudo- 
scalar, which measures the twist, linking, and writhe of mag- 
netic field lines H). A "helicity injection" into some ther- 
modynamic (or turbulent) system may create a current with 
twisting magnetic field lines. This idea has been successfully 
demonstrated in plasma experiments 11 IOU12I1 . In this case, we 
invoke a pseudo-scalar coefficient A and define a magneto- 
hydrodynamic free energy as F :— J B^ d^^x/{2flo) — XK. 
The minimizer of F gives an equilibrium magnetic field with 
a finite current (in this case, J parallels B) Ill3ll . The A 
(called Beltrami-parameter) can be interpreted as a pseudo- 
scalar chemical potential, and, introducing a grand-canonical 
ensemble of magnetic and flow fields, a Boltzmann distribu- 
tion with a finite helicity can be formulated II14I1 . 
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